Dynamics and stability of relativistic GRB blast waves

In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprecedented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov phase. Our goal is to determine the stability properties of its frontal shock. We carried out a grid-adaptive relativistic 2D hydro-simulation at extreme resolving power, following the GRB jet during the entire afterglow phase. We investigate the effect of the finite initial jet opening angle on the deceleration of the blast wave, and identify the growth of various instabilities throughout the coasting shock front. We find that during the relativistic phase, the blast wave is subject to pressure-ram pressure instabilities that ripple and fragment the frontal shock. These instabilities manifest themselves in the ultra-relativistic phase alone, remain in full agreement with causality arguments, and decay slowly to finally disappear in the near-Newtonian phase as the shell Lorentz factor drops below 3. From then on, the compression rate decreases to levels predicted to be stable by a linear analysis of the Sedov phase. Our simulations confirm previous findings that the shell also spreads laterally because a rarefaction wave slowly propagates to the jet axis, inducing a clear shell deformation from its initial spherical shape. The blast front becomes meridionally stratified, with decreasing speed from axis to jet edge.Comment: accepted for publication in A&

Transverse stability of relativistic two-component jets

Context: Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius represents a second interface where disruptions may occur. Aims: We explore the stability of stratified, rotating, relativistic two-component jets, in turn embedded in static interstellar medium. Methods: In a grid-adaptive relativistic hydrodynamic simulation with the AMRVAC code, the non-linear azimuthal stability of two-component relativistic jets is investigated. We simulate until multiple inner jet rotations have been completed. Results: We find evidence for the development of an extended shear flow layer between the two jet components, resulting from the growth of a body mode in the inner jet, Kelvin-Helmholtz surface modes at their original interface, and their nonlinear interaction. Both wave modes are excited by acoustic waves which are reflected between the symmetry axis and the interface of the two jet components. Their interaction induces the growth of near stationary, counterrotating vortices at the outer edge of the shear flow layer. The presence of a heavy external jet allows to slow down their further development, and maintain a collimated flow. At the outer jet boundary, small-scale Rayleigh-Taylor instabilities develop, without disrupting the jet configuration. Conclusion: We demonstrate that the cross-section of two-component relativistic jets, with a heavy, cold outer jet, is non-linearly stable.Comment: Accepted in A&A 24/09/200

AMRVAC and Relativistic Hydrodynamic simulations for GRB afterglow phases

We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $\gamma=100$ up to the almost Newtonian $\gamma\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.Comment: 12 pages, accepted in MNRAS, 12/01/200

Imaging a boson star at the Galactic center

Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an accretion torus around Sgr A* assuming this compact object is a boson star, i.e. an alternative to black holes within general relativity, with no event horizon and no hard surface. We show that very relativistic rotating boson stars produce images extremely similar to Kerr black holes, showing in particular shadow-like and photon-ring-like structures. This result highlights the extreme difficulty of unambiguously telling the existence of an event horizon from strong-field images.Comment: 21 pages, 9 figures, accepted in CQG; main difference wrt previous version is the last paragraph of the conclusio

Faranoff-Riley type I jet deceleration at density discontinuities "Relativistic hydrodynamics with realistic equation of state"

The deceleration mechanisms for relativistic jets in active galactic nuclei remain an open question, and in this paper we propose a model which could explain sudden jet deceleration, invoking density discontinuities. This is particularly motivated by recent indications from HYMORS. Exploiting high resolution, numerical simulations, we demonstrate that for both high and low energy jets, always at high Lorentz factor, a transition to a higher density environment can cause a significant fraction of the directed jet energy to be lost on reflection. This can explain how one-sided jet deceleration and a transition to FR I type can occur in HYMORS, which start as FR II (and remain so on the other side). For that purpose, we implemented in the relativistic hydrodynamic grid-adaptive AMRVAC code, the Synge-type equation of state introduced in the general polytropic case by Meliani et al. (2004). We present results for 10 model computations, varying the inlet Lorentz factor from 10 to 20, including uniform or decreasing density profiles, and allowing for cylindrical versus conical jet models. As long as the jet propagates through uniform media, we find that the density contrast sets most of the propagation characteristics, fully consistent with previous modeling efforts. When the jet runs into a denser medium, we find a clear distinction in the decelaration of high energy jets depending on the encountered density jump. For fairly high density contrast, the jet becomes destabilised and compressed, decelerates strongly (up to subrelativistic speeds) and can form knots. We point out differences that are found between cylindrical and conical jet models, together with dynamical details like the Richtmyer-Meshkov instabilities developing at the original contact interface.Comment: accepted in A&

Shocks in relativistic transverse stratified jets, a new paradigm for radio-loud AGN

The transverse stratification of active galactic nuclei (AGN) jets is suggested by observations and theoretical arguments, as a consequence of intrinsic properties of the central engine (accretion disc + black hole) and external medium. On the other hand, the one-component jet approaches are heavily challenged by the various observed properties of plasmoids in radio jets (knots), often associated with internal shocks. Given that such a transverse stratification plays an important role on the jets acceleration, stability, and interaction with the external medium, it should also induce internal shocks with various strengths and configurations, able to describe the observed knots behaviours. By establishing a relation between the transverse stratification of the jets, the internal shock properties, and the multiple observed AGN jet morphologies and behaviours, our aim is to provide a consistent global scheme of the various AGN jet structures. Working on a large sample of AGN radio jets monitored in very long baseline interferometry (VLBI) by the MOJAVE collaboration, we determined the consistency of a systematic association of the multiple knots with successive re-collimation shocks. We then investigated the re-collimation shock formation and the influence of different transverse stratified structures by parametrically exploring the two relativistic outflow components with the specific relativistic hydrodynamic (SRHD) code AMRVAC. We were able to link the different spectral classes of AGN with specific stratified jet characteristics, in good accordance with their VLBI radio properties and their accretion regimes.Comment: 16 pages, 12 figures, accepted for publication in A&

Circular geodesics and thick tori around rotating boson stars

Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects.Comment: Accepted for publication in CQ

No visible optical variability from a relativistic blast wave encountering a wind-termination shock

Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high performance adaptive mesh relativistic hydrodynamics code, AMRVAC, and the synchrotron emission is analyzed in detail with a separate radiation code. We find no bump in the resulting light curve, not even for very high density jumps. Furthermore, by analyzing the contributions from the different shock wave regions we are able to establish that it is essential to resolve the blast wave structure in order to make qualitatively correct predictions on the observed output and that the contribution from the reverse shock region will not stand out, even when the magnetic field is increased in this region by repeated shocks. This study resolves a controversy in recent literature.Comment: 4 figures, submitted to MNRAS letter

Two-flow magnetohydrodynamical jets around young stellar objects

We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. Our simulations show that the inner outflow is accelerated from the central object hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force. We also show that when a dense inner stellar wind occurs, the resulting disc-driven jet have a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass ejection rate in the disc-driven jet and a larger radial extension which is in better agreement with the observations besides being more consistent.Comment: Accepted for publication in Astrophysics & Space Science. Referred proceeding of the fifth Mont Stromlo Symposium Dec. 1-8 2006, Canberra, Australia. 5 pages, 3 figures. For high resolution version of the paper, please click here http://www.apc.univ-paris7.fr/~fcasse/publications.htm